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The subject of the present paper is the theoretical description of collective electronic excita- 
tions, i.e. spin waves, in the Hubbard-model. Starting with the widely used Random-Phase- 
Approximation, which combines Hartree-Fock theory with the summation of the two-particle ladder, 
we extend the theory to a more sophisticated single particle approximation, namely the Spectral- 
Density- Ansatz. Doing so we have to introduce a 'screened' Coulomb-interaction rather than the 
bare Hubbard-interaction in order to obtain physically reasonable spinwave dispersions. The dis- 
cussion following the technical procedure shows that comparison of standard RPA with our new 
approximation reduces the occurrence of a ferromagnetic phase further with respect to the phase- 
diagrams delivered by the single particle theories. 
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I. INTRODUCTION 

Recent publicationsBBiSHilBi on the theoretical de- 
scription of spin waves in the Hubbard-model demon- 
strate an increased interest in the subject. The Hubbard- 
modelfHjOijiainally introduced in 1963 by Hubbard and 
othersBt2HlJ, is one of the best investigated models in 
many-body-theory. It is commonly accepted to be the 
simplest model to hold for the itinerant as well as the lo- 
calized character of the valence electrons in partially filled 
3d-shells of transition and rare earth metals such as Ni, 
Fe, Co. Despite the continuing interest only a handful 
of exact statements can be given on the possibility of a 
ferromagnetic groundstate. 

Actually the model appears not to be very convenient 
for ferromagnetian. The exact solution of the 1 — D 
Hubbard-modelllS suggests the ground-state to be an an- 
tiferromagnetic one as well as the fact that in all dimen- 
sions the half-filled Hubbard-model can be mapped onto 
the antiferromagnetic Heisenberg-model in the strong 
coupling regime ^ > 1, where W denotes the band- 
width and U the interaction strength. Furthermore the 
Mermin- Wagner-Theoremll3 forbids the existence of mag- 
netic long range order for Z? < 3 at finite temperature T. 
The zero-bandwidth limit excludes a|farromagnetic solu- 
tion because n| — is always validtil, herein n„ is the 
occupation number of 'f(|)-spin. 

If at all, we expect a ferromagnetic phase to occur in 
the strong coupling regime > 1 , away from exact half- 
filling in I? = 3. 

The simple Hartree-Fock Approximation (HF) is able 
to deliver a ferromagnetic ground state in the Hubbard- 
model as long as the interaction strength is sufficiently 
large or respectively the density of states at the Fermi- 
surface Po{£f) is considerably high, manifested in the fa- 



mous Stoner-criterion 1 < Upo{e-p). 

It is a well known fact that HF-theory overestimates 
the existence of spontaneous ferromagnetism. In partic- 
ular the Curie-temperatures are orders of magnitude too 
high but nevertheless HF points in the right direction. 

A more sophisticated approach beside others is the 
Spectrai-IDensity-Ansatz (SDA)E3 based on Harris's and 
Lange'al3 statements who have shown that the spectral 
density is dominated by two strongly developed peaks 
at the centre of the band Tq and at To + when a fi- 
nite hopping probability is taken into account. SDA ba- 
sically is a two-pole ansatz of the form S\^a{E — fx) = 
fi'}2ii=i 2^ia-^))^i^ ~ 'E'icr(k)), whcrc the in the first 
place unknown parameters, spectral weights aia and 
peak positions i?iCT(kl_are fitted via the exactly calcula- 
ble spectral momentsEj. The SDA-approximation results 
in a much more reliable phase-diagram and more trust- 
worthy Curie-temperatures than HF-approximation. In a 
single particle theory, however, collective excitations are 
not present. The study of these excitations in itinerant 
ferromagnets requires the examination of the transver- 
sal magnetic susceptibility. Collective excitations have 
a strong influence on the stability of the ferromagnetic 
phase in particular at higher temperatures since the ex- 
citations of magnons tend to destroy the magnetic or- 
dering. This is the reason for our interest in describing 
collective excitations in the Hubbard-model. 

Facing the complexity of the many-body problem one 
usally has to apply one or more suitable and reasonable 
approximations in order to find physically interpretable 
results. Certain sequences of such approximations will be 
presented in the following sections. In Sect. |l^ we will de- 
velop the theory for the magnetic susceptibility which re- 
sults in the Random-Phase- Approximation (RPA). RPA 
combines ladder-approximation of the two-particle Green 
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function with the HF-approximation of the one-particle 
Green function and is discussed in more detail in Sect. [I]. 



In Sect. [V we will combine ladder- approximation of 



the two-particle Green function with the more sophisti- 
cated SDA-approach of the one particle Green function 
which makes a further approximation necessary, namely 
a screened interaction Vef / , but will result in physically 
reasonable behaviour of the collective mode. This is a 
new approach to the at the beginning mentioned spin 
wave problem and discussed in detail. We present our 
results in Sect. ^ and discuss them with respect to re- 
sults of different authors and will finally close with a short 
summary and an outlook in Sect. VI. 



time scale, this two-particle quantity is given by 
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{Tr [a+ (T)ak+q^ (T)a+^ (O)ap-qT (0)]) 
= xj-iEn) (5) 

The time-ordering operatoiJV in , for wich a general- 
ized Wick's theorem existsllj, sorts operators in a prod- 
uct of operators by their time-arguments. Matsubara- 
frequencies, denoted by subscript n, are given by 



II. THEORY 

A. Transversal Magnetic Susceptibility 

We use the so-called one band Hubbard-Hamiltonian 
where orbital degrees of freedom are neglected for sim- 
plicity. It reads in second quantized form in Fourier-space 
as 
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The operators aj^ (oko-) in (|l|) are the creation (annihi- 
lation) operators for fermions in k-space. 

Spin waves are collective excitations of the solid's state 
electronic system. The group of conduction electrons 
take on a coherent quantum state which is described as a 
gapless excitation with energy E = ftw(q), where q is the 
wavevector of a so-called spinwave. The spin wave dis- 
persion usually behaves as i?q ~ Dq^ for small q-values, 
where D is the stiffness constant. 

In order to describe the excited system we make use 
of the retarded polarisation propagator and will analyze 
it with help of a diagramatic technique. In our case the 
polarisation propagator is the transversal magnetic sus- 
ceptibility. It is a two-particle Green function and given 
in energy representation by 



where the spin-operators are 
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Within Matsubara formalismlli^0, which is a generalisa- 
tion of the Green function formalism to the imaginary 
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fermions 
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One obtains the retarded form of (||), the actual inter- 
esting quantity, simply by analytical continuation. The 
poles of yield the spin wave excitation spectrum. Dia- 
gramatically is represented by the graph given on the 
left side of Fig. |l|. The representation with help of the 
vertex part (shaded triangle on the right of Fig. ^) is very 
handy because it allows an approximate partial summa- 
tion. The vertex part plays a similar role as the selfenergy 



FIG. 1: Representation of the polarisation propagator with 
help of the vertex function 

in a single particle theory. The doubly lined propagators 
on the right side of Fig. |^ indicate that we have to use 
the full single particle Green function G^l^{En) instead 
of the free Green function G^'*^(ii'„) indicated by singly 
lined propagators as on the left side of Fig. ^ 

For the vertex part we have to find a practicable 
and physically reasonable approximation. Summing up 
the ladder graphs (see Fig. g) we find with help of the 








FIG. 2: Partial summation of the ladder graphs 



3 



*4 



diagram- 'dictionary 



ru(q,^„;k,0 = 1 



-1 u 



P E'' 



G^^(i?;')rn(q,i?„;p,<')- (7) 

The coupling parameter U in the Hubbard-model is a 
constant which is why we can write 

rTi(q, E„;k, E'J = T^^iq, p, <') = r(q, S„) (8) 

for the vertex part. 
Introducing 



An(q, E„) = ^ Yl G^^iK)G^+^i{En + <) 

(9) 



k E' 



we write for the whole propagator: 

^- ATi(q,^„) 



(10) 



Y+-(E ) - 1^^^' 



The abbreviation A|j(q, £"„) on the right side of (|1C 
represents the polarisation propagator of zeroth order or 
pair bubble where no interaction line connects two vertex 
points. 

Now, finding (|^) solves the problem. This task can be 
fuUfilled by first rewriting (||) with help of the spectral 
representation of the Matsubara function, i.e. we use 

oo 

G^.iEj J dE,——^^ 



The polarisation propagator then becomes A||(q, E^) = 

CO 

J f dEidE2F{E„; Ei,E2)- K^{Ei,E2). Both expression 

— oo 

( |l3| ) und (|lj) can be treated independently. 

Evaluatioa of (|l|) requires the sum over Matsubara 
frequencies^ which results in 



F{En;Ei,E2) = 



f.{E2)- f.{E^) 
iEn - {E2 - El) 



(15) 



where /-(^^i.s) = ^mhr^- 

The polarisation propagator A|| (q, i?„) now looks with 
Ei,2 — > (£'i,2 - /i) like 
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Here /_(£^i_2) in the numerator is the Fermi function. 

In order to represent the transversal susceptibility (|^) 
as a functional of a product of two single particle spec- 
tral densities we performed a ladder approximation which 
delivered the desired form (|l^). Up to now we did not 
make any statement about the one-particle spectral den- 
sities in (^6|) . On account of the complexity of the many- 
body problem we have to apply a suitable approximation 
for the single particle spectral densities in (|l6|). In the 
following section we will do exactly this with the HF- 
approximation. 



Gk+qliK+En) - J dE2-^ 
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and write for A|| (q, £"„ 
-1 1 
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Au(q,£^n) = ^^EE / J dE2X 



Skf{El) S'k+qx(i?2) 



iE^-Ei i{E'^ + E^)~E2 
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For technical reasons we introduce two other abbrevia- 
tions namely. 



i{E',^+En)- E2 



and 



Kq{Ei,E2) = ^E^kT(i?i)W(^2). (14) 



III. RANDOM-PHASE-APPROXIMATION 

Evaluation of (|l^) requires a three-dimensional sum- 
mation over wavevectors k and k -I- q respectively. 

The energy-dispersion for the conduction band in 
Tight-Binding- Approximation is given for the 6cc-lattice 
by 

g^j^^&cc _ ZiTicos{kxa)cos{kya)cos{kza) + Tq (17) 
and for the /cc-lattice by 

e(k)/^^ = AZiTi[cosih,a)cos{hya) 



+cos(-kya)cos(-kza) 
2 2 



+cos{^kza)cos{^k^a)]+To. (18) 

(13) The Hopping-integral Ti for next nearest neighbours is 
given by Ti = — , where we used a bandwidth of = 
leV for all numerical calculations, Zi denotes the number 
of next nearest neig hbours {Z^""^ = 8 and Zl"" = 12), the 
center of the Band was choosen to be Tq = 0. In order 
to reduce the numerical effort to a minimum we use a 
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similar q-separation procedure as given in0 where the 
remaining k-summation can be transformed into a one- 
dimensional energy-integration. This is due to the fact 
that the expression 



lF(e(k)) = I F{x)p 



BDOS 



where 



(19) 



(20) 



is true for any k-dependent function whose k- 
dependence is only given by a e(k)-dependcnce. 

This procedure alloms us to apply the numerical for- 
mulas given by Jelittoll3 for the Bloch density of states 
(BDOS). Picture I shows the BDOS for the fee and bcc- 
lattice resulting from these formulas. 
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FIG. 3: Bloch Density of States, bandwidth W — leV, center 
of the band To = 

For the situation where q = 0, however, the q- 
separation is not needed and one easily can check if there 
exsists a gapless excitation with cj(q) — > as proposed 
by LangeEil. Useful is the following representation of the 
HF-spectral density 

SkaiE)^hSiE-eik)-U{n^,)+^i). (21) 



The propagator (g^ reduces with (|2l| ) and i?k,cr = Ck 
U{n„) to 



A 



n 



/_(£;k+q,T)-/-(i?k,i) 



(22) 



For q = one easily can perform the summation over 
wave vectors because 



^E/-(^k,.) = ^E(-'^(k))- (23) 



The condition for the denominator of ( p^ ) to be zero is 
l-C/-ReAti(q = 0,£;) = (24) 



or equivalently 



1 



-U ■ 2m 



E-U ■2m + iO+ 







(25) 



Here m — ^{ni — ni) denotes the magnetisation. One 
immediatly sees that this condition is fuUfilled only for 
E = which means that using HF-theory a well defined 
gapless excitation exsists. 

EdwardscB already stated that excited states with one 
reversed spin in bandferromagnets must have an energy 
of the form E ^ e(k') — e(k) + J where J is a k- 
independent constant that contains the correlation en- 
ergy. Clearly in HF-theory this is the case. In the 
event of transitions between majority- and minority-spin 
band the excitation energy is given by A£^||(k;q) = 
e(k -I- q) — e(k) -I- AE^x, where AEex = 2mU is the rigid 
exchange splitting energy of the HF-thcory. 



IV. MODIFIED RANDOM PHASE 
APPROXIMATION 

In this section we will use the SDA-theory as the 
single particle input for the two particle propagator 
(p6|). Compared to HF- approximation SDA is a real 
many-body theory able to deliver ferromagnetism in 
the Hubbard-model and tested by means of compar- 
ison to different advanced many-body theories con- 
cerning ferromagnetism in the Hubbard-model. Com- 
parison of SDA, Modified-AUoy-Analogy (MAA) and 
Modified-Perturbation-Theorie with the numerical ex- 
act Quantum-Monte-Carlo-Simulationc^ indicates that 
the qualitative behaviour of the SDA phase-diagram 
is correctly displayed. Furthermore are the Curie- 
temperatures orders of magnitude lower than in HF- 
theory and therefore much more reliable. HF-theory 
overestimates the occurence of ferromagnetism in the 
Hubbard-model, although the tendency is right, in a in- 
credulous manner. 

Though one cannot simply check the validity of the 
theory by putting the SDA-selfenergy (or equivalently 
the spectral density) for q = into the expression ( p^ ) 
because the selfenergy in SDA consists of a Stoner or 
HF-like part and a spin dependent rest. This spin- 
dependence prevents us from being able to check the va- 
lidity in a similar manner as with HF-theory. 



^SDA 
■"kCT 



(E) = Un^, + 



E- B- 



(26) 



The spin dependence of the second term of (^6|) is mainly 
given by the so called band shift B-^- This spin- 
dependence in SDA-theory is the main improvement and 
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leads to a shift of and |-quasipartical-density of states 
(QDOS), for further details on the SDA se£i. 

However, the similarity of ( p6| ) to the Stoner self energy 
let us hope that connecting SDA with ladder approxima- 
tion could work out and might even hold an improvement. 
But numerical evaluation of ( [l6| ) using SDA-theory as 
input in form of a sum of weighted (5-functions shows 
negative excitation energies which are unphysical. 

With no apparent reason both approximations, the 
ladder approximation for the two-particle propagator and 
the SDA for the single particle propagator are not fully 
compatible. 

The idea now was to introduce an effective Coulomb 
interaction Vefj in order to replace the Hubbard-U in 
( |lo| ) and hopefully be able to shift the negative excitation 
energies to the right energetic position. This was done 
numerically via the condition 



Veff = 



u 



Re 



An(O,0) 



i/ReA|^(q = 0,^ = 0) = ^ 
else 



Indeed, replacing the Hubbard-C/ in the polarisation 
propagator ( ]lO| ) with Veff given by the above condition 
the unphysical behaviour could be repaired. The suscep- 
tibility (10) now is to look at as 



1- 14// •A;^^^^(q, £;„,£/) 



(27) 



The term A^"^^^) in (|^) is with the appropriate 
input of the SDA-single particle density of states. 

The introduction of the effective or screened interac- 
tion of course is to rate as yet another approximation. 
But the interaction in a two particle theory not neces- 
sarily has to be the same as the interaction in a single 
particle theory. It is the essence of a two particle theory 
to take additional correlations into account. HF-theory 
for instance neglects the correlation between particles 
of opposite spin. This can be seen by the diagramati- 
cal representation of the self-energy in HF-theory, where 
the exchange-interaction part does not contribute to the 
self-energy. The only contribution to the HF-self-energy 
comes from the direct Coulomb-term, wich describes the 
interaction of an electron with an effective field caused by 
the other particles in the system. The contribution com- 
ing from the exchange-interaction, however, is taken into 
account in ladder approximation as can be seen from (|^). 
Diagramatically spoken we did nothing else than another 
renormalization in order to be able to use SDA-theory. 
The broken lines in Fig. ^ now has to be imagined as 
broken double lines indicating an effective interaction in 
the vertex. 

The effective potential as a function of the bare 
Hubbard-C/ can be seen in Fig. || for the /cc-lattice at 
T — OK and three different band-occupations. Since the 
effective potential seems to saturate for larger values of 
U we restricted ourselves to f/-values up to 7eV. 

The SDA-QDOS is mainly given by two wheighted S- 
like peaks. This two-peak structure is conserved in our 




FIG. 4: Veff ~ f{U) for band occupations n = 1.4, n — 1.6, 
n = 1.8 at T = OK 



MRPA-method. This can be seen from Fig.|^, where 

Rext~^^"^''^\En) and Re{l-Vef fAui^l, En)) are plot- 
ted for q = and q = vr as a function of Energy. The 
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FIG. 5; Upper graphs: Rex^~'''^' ^^^\En) for q = and q = 
(7r,7r,7r), the two-peak structure of the two-particle spectral 
density is clearly visible in the upper left graph; the lower 
grap h shows Re{l — Ve//A||(q, iS^)) to indicate the poles of 



effect of the SDA is that the upper set of q-dependent 
curves, i.e. the one shifted by '--^ C/ as aresult of the SDA- 
input, shows a dispersion too. This set could be identified 
as optical modes which are not preseut, in RPAj.heory. 
Calculations for the real substances FecJ and Nic3, how- 
ever, indicate that a separation between acoustical and 
optical modes about U is one order of magnitude to 
large. That is why we only did look at the low lying 
excitation mode. 
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V. RESULTS AND DISCUSSION 

In this section we present the resuhs of MRPA and 
RPA. The latter one we used as the estabhshed 'test' 
theory for our purposes. 

In MRPA we restricted ourselfs to the /cc-lattice since 
an asymmetric quasi particle density of states (Q-DOS) 
supports the tendency towards a ferroiBagnetic aligne- 
ment of the spins in the Hubbard-modelE^. 

Figure ^ shows isotropic spin wave dispersions as a 
function of U for the /cc-lattice in MRPA. The disper- 
sion curves were computed along different directions of 
high symmetry in the First-Brillouin-Zone (FBZ). This 
isotropy for the /cc-lattice is in jagreement with results 
observed using Gutzwiller-ansatzEJ. 




«.RPAn=0.7bcc(100) 
™ MRPA n=1.3 fee (100) 




2 u [eV] 4 6 

FIG. 7: Behaviour of magnon energies at the edge of the FBZ 
as a function of U for RPA and MRPA at T = OK, qualitativ 
similar behaviour for larger [/-values in both theories 



diagram which is to see in Fig. ^ 

Figure || compares HF-phase diagrams with the corre- 
sponding phases where the spin wave mode, indicated by 
subscript SW, exsists for three interaction strengths U for 
the 6cc-lattice. One example shows a similar behaviour 
for the /cc-lattice. 

Numerically we distinguished the SW phase by looking 
at the stiffness or the magnon energy at the edge of the 
FBZ respectively. This means that whenever the stiffness 
D becomes negative the ferromagnetic state breaks down. 



FIG. 6: Magnon dispersion in MRPA along different high 
symmetry directions of the FBZ for the /cc-lattice, param- 
eter: U — 2.71ey ; 5.00ey ; 6.50ey at fixed band occupation 
n = 1.4, T — OK, see the bordered box for values of Ve// 

Figure || also shows the behaviour of the stiffness D 
in dependency of the correlation parameter U. For in- 
creasing values of U the magnon energy at the edge of 
the FBZ weakens or in other words D decreases with 
increasing U. In Fig. £'(q = (tTjTT, tt)) = f{U) is plot- 
ted for RPA and MRPA for fixed band occupations in 
either case. We see that in RPA theory the magnon en- 
ergy at the edge of the FBZ continuously increases up 
to a maximum and afterwards decreases with higher val- 
ues of U. In MRPA the magnon energy at the edge of 
the FBZ decreases continously with increasing U. Thig 
is in contrast to a statement given by Biinemann et.al.B 
where the authors found the opposite behaviour using 
(multiband) RPA namely increasing D with increasing 
U. The behaviour of decreasing D with increasing U on 
the other hand was observed by Biinemann et.aljj using 
their recently developed multiband Gutzwiller-ansatz. 

We rather think that the decreasing D with increasing 
?7-behaviour is an intrinsic feature of the Hubbard-model. 

To us the most surprising result is that using RPA we 
could compute a qualitatively rather new Tc{n)-p]iase 
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FIG. 8: RPA: Tc(n)-phase diagrams, the upper curves are 
Tc (?i)-phase diagrams computed with single-particle HF- 
theory denoted by subscript HF, in each case the lower curves 
indicate the phases where spin waves are excited denoted by 
subscript SW 

The dramatic influence of the spin wave mode on 
the stability of the ferromagnetic phase is obvious from 
Fig. ||. Since Stoner-theory overestimates the possibility 
of ferromagnetic ordering in the Hubbard-model the lack 
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of the single particle theory is compensated by the cor- 
responding two-particle theory and even results in new 
Curie-temperatures. 

This feature of thej-|RPA is also described for the case 
of Neel-temperaturea3. 

In the modified theory MRPA the influence of the col- 
lective mode on the shape of the rc(r7,)-phase diagram 
can be seen from Fig. |9|. The resulting phase diagram 
again is physically understandable. The excitation of 




FIG. 9: MRPA: Comparison of single and two-particle phase 
diagram, in each case the upper curve is the rc(n)-curve from 
SDA-theory denoted by subscript SDA, the lower curves are 
the 'collective' phases denoted by subscript SW, fcc-lattice 

spin waves reduces Curie-temperatures and the region of 
band occupation but the qualitative shape more or less 
stays the same as in the single particle theory. 

The ferromagnetic phase for both theories in general 
is reduced due to weakening of the spin wave disper- 
sion close to half filling and the temperatur driven in- 
creased magnon excitation. From figure ^ can be seen 
that the Tc-curves have their peaks for n 1.4. This 
feature is due to the fact that the bandoccupation that 
favours a ferromagnetic alignment of the spins in the 
SDA-approach is mainly given by the lattice-structure. 
For the fcc-lattice the favourable bandoccupation is n > 1 
and for n ^ 1.4 the ferromagnetic groundstate is almost 
saturized. The concrete value for Tq is mainly given as 
a function of the correlation strength U . 



Comparison with Tc(n)-ph|ase diagrams of different 
single particle approximationsES including SDA and tha 
numerically exact Quantum-Monte-Carlo-SimulationllB 
demonstrate that besides Curie-temperatures the region 
of band occupation where ferromagnetism exists is re- 
duced. This indicates that our 'coUective'-phase dia- 
grams are much more reliable then the ones computed 
using the corresponding single particle theories. 

VI. SUMMARY AND OUTLOOK 

In the present paper we computed spin wave disper- 
sions and Tc (n)-phase diagrams using standard RPA the- 
ory and a modified theory where a reliable approximation 
for a single particle Green function is used as an input 
for the two particle Green function. Although we did not 
follow a stringent mathematical procedure for pur MRPA 
as described exemplary for the case of RPA incJ we were 
able to show that our results using the modified theory 
MRPA gives physically reasonable results. These results 
are comparable to results which we get using standard 
RPA and to results of a different method namely the 
Gutzwiller-ansatz as described|-|iiM,-|rAlthough RPA is a 
widely used standard techniquea'BoBQlD to our knowledge 
a 'complete' Tc(n)-phase diagram was not published be- 
fore. 

Our aim was to get a deeper insight into the phe- 
nomenon of collective excitations in a model bandfer- 
romagnct. To a 'first' approximation our results show 
that just counting on a single particle theory in order to 
describe the problem of collective magnetism is not infor- 
mative enough. Collective excitations strongly influence 
the stability of a ferromagnetic state and one might use 
'collective'-phase diagrams as tests for single particle the- 
ories. 

If one wanted to improve the whole theory one had 
to find an approximation for the vertex function bet- 
ter, i.e. more informative, than ladder approximation 
and an approximation for the single particle selfenergy 
at the same theoretical level. Parquet-Ansatz in con- 
nection with Coherent-Potential- Approximation CPA as 
described for the electrical conductivity in the Anderson 
modelnJ might be a good starting point. A CPA- type 
approximation be able to hold for ferromagnetism in this 
respect is the Modified-AUoy- Analogy (MAA)cJ which 
already yields an iaaproved phase-diagram compared to 
the SDA-approachEl. 
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